Structure Function Scaling in Compressible Super-Alfvenic MHD Turbulence 
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Supersonic turbulent flows of magnetized gas are believed to play an important role in the dynam- 
ics of star-forming clouds in galaxies. Understanding statistical properties of such flows is crucial 
for developing a theory of star formation. In this letter we propose a unified approach for obtaining 
the velocity scaling in compressible and super-Alfvenic turbulence, valid for arbitrary sonic Mach 
number, Ms. We demonstrate with numerical simulations that the scaling can be described with 
the She-Leveque formalism, where only one parameter, interpreted as the Hausdorff dimension of 
the most intense dissipative structures, needs to be varied as a function of Ms. Our results thus 
provide a method for obtaining the velocity scaling in interstellar clouds once their Mach numbers 
have been inferred from observations. 

PACS numbers: 98.38.Am, 95.30.Qd, 98.38.Dq, 95.30.Lz 
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1. Introduction. Star formation occurs in our galaxy 
primarily inside cold interstellar clouds. Observations in- 
dicate that gas motion inside these clouds is highly super- 
sonic in the broad range of scales 0.1 — 1000 pc. At these 
scales the clouds lack any characteristic structure and 
their scaling properties are to some extent universal. This 
universality is the signature of the supersonic and super- 
Alfvenic magneto-hydrodynamic (MHD) turbulence that 
governs the dynamics of the clouds Q, 0, 0, 0, H, IE ■ 

The properties of turbulent MHD flows are controlled 
by the rms sonic Mach number, Ms, (ratio of rms flow 
velocity and sound velocity) and the rms Alfvenic Mach 
number, Ma, (ratio of rms flow velocity and Alfven ve- 
locity). In interstellar clouds these numbers are large, 
Mg ~ 1 — 30, Ma > 1. The corresponding supersonic 
motionproduces the observed strong density enhance- 
ments p , some of which are Jeans-unstable and undergo 
gravitational collapse. This initial 'turbulent' stage of 
density fragmentation should be described statistically, 
in terms of velocity and density probability distribution 
functions or correlators; its properties are inherently re- 
lated to the statistical properties of compressible MHD 
turbulence Q. 

Early studies suggesting the importance of astro- 
physical supersonic turbulence were undertaken by von 
Weizsacker and later by Larson who discovered 
that the spectrum of the velocity field in molecular clouds 
was steeper than the spectrum predicted by the standard 
Kolmogorov model. The discrepancy was attributed to 
the strong intermittency of supersonic turbulence. How- 
ever, the analytical understanding of such turbulence re- 
mained rather poor. 

Advances in computational and observational tech- 
niques have recently allowed measurements of power 
spectra and higher order correlators of velocity and den- 
sity fields, reviving the interest in supersonic turbulence. 
Modern observations generate surveys covering a range 



of scales comparable or even exceeding the resolution of 
numerical simulations. The interstellar medium can now 
be viewed as a laboratory for investigating supersonic 
motion unachievable with Earth-based experiments. 

The turbulent field can be characterized by the scaling 
of its structure functions [see below] . The intermittency 
is defined as the departure of this scaling from the Kol- 
mogorov one. So far, only the limiting cases of Ms 3> 1 
and Ms = have been investigated. An analytical de- 
scription of supersonic turbulence, Ms 3> 1, has been 
proposed in |l2( and confirmed numerically in |13lll4Lll5| . 
It is based on the observation that the most intense dissi- 
pative structures in supersonic flows are two-dimensional 
shocks. This approach makes use of the model by She 
and Leveque |l6| for incompressible turbulence, Ms = 0, 
where the dissipative structures are one-dimensional vor- 
tex filaments. 

In this letter we report the results of a series of numer- 
ical MHD simulations where we traced the change of the 
velocity structure functions over the broad range of Mach 
numbers 0.3 < Ms < 10. We discovered that their scal- 
ings can be described by the general model where only 
one parameter, the dimension of the most intense dis- 
sipative structures, needs to be changed with Mg. The 
scalings obtained in 0, ^| naturally appear as the lim- 
iting cases of this unifying model. Our results suggest 
that the velocity scaling within interstellar clouds can be 
inferred from measurements of their Mach number. Ob- 
servationally, this is a much easier task than measuring 
the structure function directly. We present the main re- 
sults here; the detailed discussion will appear elsewhere. 

2. Velocity Structure Function Scaling. The velocity 
structure functions of order p are defined as: 



S p (£) = {\u(x + £) - u(x)\ p ) oc^ p ) 



(1) 



where the velocity component u is parallel (longitudinal 
structure function) or perpendicular (transversal struc- 
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FIG. 1: Compensated velocity structure functions of order 
p = 1 to p — 10 (bottom to top) from the Ms=9.5 experiment. 
To avoid confusion, the structure functions of order p = 1 and 
p — 2 have been multiplied by 0.2 and 0.7 respectively and 
only the upper part of error bars larger than 100% is shown. 



FIG. 2: Structure function exponents normalized to the values 
predicted by |T^| for supersonic turbulence (dashed-dotted 
line). Values from five experiments are shown, from Ms=9.5 
(bottom) to Ms=0.4 (top). The Ms=0.8 case is not shown 
as it almost overlaps with the Ms =0.4 experiment. 



ture function) to the vector £ and the spatial average is 
over all values of the position x. 

The scaling of velocity structure functions in incom- 
pressible turbulence is best described by the She-Leveque 
formula [T^ . fl7 | : 



C(p)/C(3)=p/9 + 2 l-(2/3) 



\P/3 



(2) 



This is a specific case of the more general formula 

Cb)/C(3)-7P + C(l-/3 p ) (3) 



with the normalization condition, 

C = (l-3 7 )/(l-/3 3 ). 



(4) 



Two basic assumptions are used to derive (0, as dis- 
cussed in [l8| . The first assumption is the existence of a 
symmetric relation between the intensity of fluctuations 
F p (£) and F p+1 {£) with respect to a translation in p: 



F p+1 (£) = A p F p {£fF {oo) {£) 1 



(5) 



where the pth order intensity of fluctuations is defined 
as Fp(£) — S p+ i(£)/S p (£), and A p are constants indepen- 
dent of £ (generally found to be independent of p as well) . 
The functions F p (£) are associated with higher intensity 
fluctuations with increasing values of p. The parameter 
(3 is a measure of intermittency. In Kolmogorov's 1941 
model [l^, P ~* 1 an d Cip) = p/3 This is the limit of no 
intermittency. Lower (3 corresponds to higher degree of 
intermittency, or to the increasing role played by strong 
persistent fluctuations (structures). 



The second assumption is 



F n 



07 
<->3 • 



(6) 



The value of the parameter 7 is related to very high or- 
der moments and is therefore difficult to constrain with 
experiments. Kolmogorov's turbulence corresponds to 
7=1/3. 

The parameter C has been interpreted as the Haus- 
dorff codimension of the support of the most singular 
dissipative structures |l7|. In incompressible turbulence 
the most dissipative structures are along vortex fila- 
ments with Hausdorff dimension D = 1 and so C = 2. 
Boldyrev |3| has proposed to extend the scaling © 
to supersonic turbulence, with the assumption that the 
Hausdorff dimension of the support of the most singular 
dissipative structures is D = 2 (C = 1), because dissipa- 
tion of supersonic turbulence occurs mainly in sheet-like 
shocks. He used 7 = 1/9 as in incompressible turbulence, 
which implies /3 3 = 1/3, that is a higher degree of inter- 
mittency than in incompressible turbulence. With these 
parameters one obtains the velocity scaling proposed for 
supersonic turbulence in |l2|: 



C(p)/C(3)=p/9 + l-(l/3) 



P/3 



(7) 



This velocity scaling has been found to provide a very 
accurate prediction for numerical simulations of highly 
supersonic and super- Alfvenic turbulence [Hi . Il^j . The 
same scaling has also been proposed for incompressible 
MHD turbulence |2fj| , where dissipation occurs mainly in 
two dimensional current sheets. 
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FIG. 3: "f3 test" plot for the models with Ms=0.4 and FIG. 4: Hausdorff dimension of the support of the most dis- 



M s =9.5. The M s =0.4 model is shifted upward by 0.2 y- 
units. The solid lines are least square fits with slope j3. 



sipative structures versus rms sonic Mach number. 



In the present work we obtain the structure func- 
tion scaling in simulations of compressible super- Alfvenic 
MHD turbulence as a function of the rms sonic Mach 
number of the flow, Ms, for constant value of the rms 
Alfvenic Mach number, Ma> 1 (the limiting case of 
strong magnetic field was investigated analytically and 
numerically in 0, H2, |23| ) . We use an isothermal equa- 
tion of state and vary the value of Mg in different exper- 
iments by varying the thermal energy. As a result, the 
initial value of Ma=10 remains unchanged from run to 
run. In our numerical method the total magnetic flux is 
conserved. However, the magnetic energy is amplified. 

The initial value of the ratio of average magnetic and 
dynamic pressures is (P m )in/(Pd)in = 0.005 for all runs. 

The value of the same ratio, averaged over the time 
interval for which the structure functions are computed, 
is given in Table 1. Because the average dynamic pres- 
sure is always in excess of the average magnetic pres- 
sure, (P m )/(Pd) ^ 1) the dissipation is expected to occur 
primarily in vortex filaments (roughly corresponding to 
magnetic flux tubes) in subsonic runs and in shocks in 
supersonic runs. In the next section we argue that when 
we increase the value of Ms, the Hausdorff dimension 
of the support of the most dissipative structures grows 
continuously from D«l to D » 2. 

3. Results. We have solved the three dimensional 
compressible MHD equations, for an isothermal gas with 
initially uniform density and magnetic fields, in a stag- 
gered mesh with 250 3 computational cells and periodic 
boundary conditions. Turbulence is set up as an initial 
large scale random and solenoidal velocity field (gener- 
ated in Fourier space with power only in the range of 
wavenumber 1 < K < 2) and maintained with an exter- 



nal large scale random and solenoidal force, correlated 
at the largest scale turn-over time. Details about the 
numerical method are given in Q. 

We have run six experiments with different values of 
Ms, in the range 0.4 < M s < 9.5 (Table HJ), for approx- 
imately 10 dynamical times. The transverse structure 
functions of velocity are computed for the last four dy- 
namical times of each experiment, up to the order p = 10. 
Their power law slope, C(p)i i s obtained from a least 
square fit in the range 2 < £ < 96. The standard devia- 
tion of the slope is used as an estimate of the uncertainty 
ofC(p). 

Figure ^ shows the structure functions plotted versus 
the third order one for the Ms=9.5 experiment. The 
structure functions are compensated by their least square 
fit slope, ((p). We study the structure functions rela- 
tive to the third order to exploit the concept of extended 
self-similarity 0, 0] , that is the property of the rela- 
tive scaling to extend up to a range of scales affected by 
dissipation (numerical dissipation in this case). 

The exponents normalized to the third order one, 
are plotted in Figure |2 divided by their value predicted 
in :i2| and given by Q. The M s =0.4 and M s =0.8 mod- 
els follow the She-Leveque scaling of incompressible tur- 
bulence. As the value of Ms increases, the velocity scal- 
ing becomes gradually more intermittent. At Ms=9.5 
the scaling is very close to the analytical prediction for 
supersonic turbulence , as previously shown 0, 0] . 

The validity of the assumption of hierarchical symme- 
try expressed by (J5J, can be tested with a log-log plot 
of F p+ i/F2 versus F p /Fi, obtained by varying the argu- 
ment £ in If the plot is a straight line, the parameters 
A p are constants independent of £ and the assumed hi- 
erarchical symmetry expressed by (jSJ) is satisfied. This 
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(Ms) 


0.39 


0.77 


1.5 


2.6 


4.8 


9.5 


(P m )i„/(Pd)in 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


{Pm)/(Pd) 


0.15 


0.14 


0.09 


0.11 


0.09 


0.08 


P 


0.87 


0.87 


0.84 


0.79 


0.73 


0.69 


D 


1.0 


1.0 


1.34 


1.70 


1.90 


2.01 



TABLE I: Values of the rms sonic Mach number, average mag- 
netic to dynamic pressure ratio at initial time and averaged 
over the time interval for which the structure functions are 
computed, f3 parameter from fitting C(p)/C(3) and Hausdorff 
dimension of the support of the most dissipative structures. 
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is called the 11 (3 test" in [18J. Our numerical simulations 
pass the (3 test for all values of Ms. In Figure |3| we show 
the two extreme cases, Ms=0.4 and Ms=9.5. We have 
used values of F p (£) for all values of i used in the least 
square fit estimate of the exponents C(p). I n ai l the runs 
we find that the constants A p in J^J are independent not 
only of £, but also of p, as in previous numerical and lab- 
oratory studies 0, |2f| H|| . The plot in Figure |3] is in 
fact obtained varying both t and p. 

The slope of the plot in Figure provides an estimate 
of j3. For the models with Ms from 0.4 to 9.5 we obtain 
(3 = 0.82, 0.81, 0.75, 0.74, 0.77 and 0.77. We instead find 
no evidence of a systematic variation in the parameter 7 
with Ms and our data is consistent with 7 = 1/9. 

The (3 test gives support to the basic assumption of the 
She-Leveque formalism, but does not necessarily yield 
the value of (3 providing the best fit to the structure 
function scaling ((p)/((3). Values of (3 computed from a 
best fit to C(p)/C(3) from our numerical experiments are 
slightly different from the (3 test values (see Table [|J. We 
use these best-fit (3 values to compute the parameter D. 

^From the definition of the parameter C in (@J), the 
values of (3 estimated above and 7 = 1/9 we can compute 
D = 3 — C as a function of Ms. The values of (3 and D 
are given m Table [Q The value of D grows with Ms, 
from D w 1.0 in the subsonic runs to D w 2.0 in the 
most supersonic runs. 

4. Conclusions. We propose a unifying description of 
the velocity structure function scaling of compressible 
and super-Alfvenic MHD turbulence for an arbitrary 
value of the rms sonic Mach number of the flow, Ms. 
Only one model parameter, interpreted as the Hausdorff 
dimension of the support of the most dissipative struc- 
tures in the She-Leveque formalism, needs to be varied 
as a function of Ms. Besides its general importance, 
our result implies that the velocity structure function 
scaling in interstellar clouds can be predicted because 
the value of Ms in the clouds is easily estimated from 
observations. This may help improve our understanding 
of the process of cloud fragmentation that leads to the 
birth of stars. 
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